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ABSTRACT 

During an extensive multiwavelength campaign that we performed in 2013-14 the prototypical Seyfert 1 galaxy NGC 5548 has been 
found in an unusual condition of heavy and persistent obscuration. The newly discovered “obscurer” absorbs most of the soft X- 
ray continuum along our line of sight and lowers the ionizing luminosity received by the classical warm absorber. Here we present 
the analysis of the high resolution X-ray spectra collected with XMM-Newton and Chandra throughout the campaign, which are 
suitable to investigate the variability of both the obscurer and the classical warm absorber. The time separation between these X-ray 
observations range from 2 days to 8 months. On these timescales the obscurer is variable both in column density and in covering 
fraction. This is consistent with the picture of a patchy wind. The most significant variation occurred in September 2013 when the 
source brightened for two weeks. A higher and steeper intrinsic continuum and a lower obscurer covering fraction are both required 
to explain the spectral shape during the flare. We suggest that a geometrical change of the soft X-ray source behind the obscurer cause 
the observed drop in the covering fraction. Due to the higher soft X-ray continuum level the September 2013 Chandra spectrum is 
the only X ray spectrum of the campaign where individual features of the warm absorber could be detected. The spectrum shows 
absorption from Fe-UTA, O iv, and O v, consistent to belong to the lower-ionization counterpart of the historical NGC 5548 warm 
absorber. Hence, we confirm that the warm absorber has responded to the drop in the ionizing luminosity caused by the obscurer. 

Key words, galaxies: individual NGC 5548 - galaxies: absorption lines - X-rays: galaxies 


1. Introduction 

In the X-ray band, active galactic nuclei (AGN) are variable 
emitters. The origin of this variability, that can be even large 
and fast (e.g., Matt et al. 2003), is not fully understood yet 
(see e.g. McHardy et al. 2006; Turner et al. 2009; Ponti et al. 
2012). Variable absorption of the X-ray radiation along the 
line of sight is one of the possible explanations. Indeed, many 
absorbing components, spanning a broad range in ionization, 
can be detected in AGN X-ray spectra. 

Cold neutral absorption is able to strongly suppress the soft 
X-ray flux and to change the curvature of the X-ray spectrum. 
In recent years, evidence for variability due to cold X-ray 


absorption in both type 1 and type 2 AGN has increased. In 
many cases changes in the absorber column density and/or 
covering fraction on a few hours-few years timescale have 
been reported e.g., in NGC 4388 (Elvis et al. 2004), NGC 4151 
(Puccetti et al. 2007), and 1H 0557-385 (Coffey et al. 2014). On 
these same timescales, the X-ray absorbing column density may 
even drastically switch from Compton thick to Compton thin 
as observed in several type 2 AGN e.g. NGC 7582 (Piconcelli 
et al. 2007), UGC 4203 (Risaliti et al. 2010), and NGC 454 
(Marchese et al. 2012). This rich phenomenology suggests that 
cold gas is present even in the innermost region of AGN (see 
Bianchi et al. 2012). This material is probably patchy and may 
belong to the Broad Line Region (BLR) or to a clumpy torus 
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Fig. 1 . From the top to the bottom panel: the observed light curves 
of NGC 5548 in two X-ray bands and in the UV (4=2030 A). These 
curves are obtained from the daily Swift monitoring performed during 
our campaign. In each panel, the pink horizontal dashed line marks the 
flux level measured by Swift at unobscured epochs (2005 and 2007). 
From the left to the right, the vertical lines indicate the first and the last 
'XMM-Newton observation of summer 2013 (dashed lines), the three 
Chandra observations of September 2013 (dotted lines), and the last 
two observations of the XMM -Newton program (in December 2013 
and February 2014, dashed lines). 


(e.g., Miniutti et al. 2014). In the best studied Seyfert 1.8 source 
NGC 1365 a single cloud has been monitored while eclipsing 
the central X-ray source for a few hours (Risaliti et al. 2009). In 
the last decade, this source has been observed in a vast variety 
of spectral appearances, ranging from Compton thick (Risaliti 
et al. 2005) to an almost unobscured condition more reminiscent 
of a pure type 1 AGN (Walton et al. 2014). 

Absorbers at higher ionization (the so-called warm absorbers, 
WA) imprint discrete absorption lines on ~50% of type 1 
AGN spectra (Crenshaw et al. 2003b). These features, that fall 
mainly in the soft X-ray (e.g. from Nen-Nex, Fex-Fexxiv, 
and Oiv-Ovii) and in the UV (e.g., Ovi, Mgn, Cii-Civ, and 
Si n—Si iv) domain, are usually blueshifted with respect to the 
systemic velocity, indicating therefore a global outflow of the 
absorbing gas (Crenshaw et al. 2003b). In local Seyfert 1 galax¬ 
ies absorption lines have usually a narrow profile (FWHM~ 10 2 
km s' 1 ). In about 15% of optically selected quasars (Weymann 
et al. 1991) broad absorption lines (BAL) with a typical width 
of 10 4 km s' 1 have been also observed (Hamann & Sabra 
2004). In the last ten years observational campaigns providing 
simultaneous high resolution UV and X-ray spectroscopy have 
been performed for a handful of local AGN (see Costantini 
2010, for a review). In some cases (e.g., NGC 4151, Kraemer 
et al. (2005) Mrk 279, Arav et al. (2007); Costantini et al. (2007) 
1H 0419-577 Di Gesu et al. (2013)), it has been established 
that X-ray and UV narrow absorption lines (NAL) may be the 
manifestation of the same outflow, even though it is not trivial to 
determine how the X-ray and UV absorbing gas are physically 
and geometrically related. The UV and X-ray absorbing gas 
could be co-located (e.g., Mrk 509, Ebrero et al. 2011), with 
the UV absorbing components being possibly denser clumps 
embedded in a more highly ionized wind (e.g., NGC 4051, 
Kraemer et al. 2012). WAs are usually complex multicomponent 
winds spanning a broad range in velocity and in ionization. 


The lower ionization components produce the UV lines while 
the higher ionization phases are seen only in the X-rays. An 
intermediate phase producing absorption lines in both bands 
may in some cases be present. 

Variability in the warm absorption may in principle contribute 
to the overall X-ray variability of AGN on different timescales. 
In order to assess the WA variability on long (e.g., ~years) 
timescales it is necessary that high quality multiepoch spec¬ 
troscopy is available, which is seldom the case. A multiepoch 
study of the WA has been attempted for instance in Mrk 279 
(Ebrero et al. 2010) without finding significant variability. In 
a different case (Mrk 841, Longinotti et al. 2010), comparing 
two different observations taken ~4 years apart, a moderate 
decrease in the WA ionization as the continuum dims has been 
observed. An even more noticeable long-term WA variability 
has been reported in the case of Mrk 335 (Longinotti et al. 
2013), where the emergence of an ionized outflow that was 
not historically present has been observed in 2009. On shorter 
timescales, changes in the WA opacity or ionization have been 
observed for instance in NGC 3783 (~31 days, Krongold et al. 
2005) and NGC 4051 (few ks-few months, Krongold et al. 
2007; Steenbrugge et al. 2009). 

The timescale over which absorption lines are observed to vary 
can be used to measure the distance of the absorbing gas from 
the ionizing source (see Crenshaw et al. 2003b). Indeed, for 
photoionized gas in equilibrium, the recombination timescale 
depends on the gas number density n. Besides the distance r it is 
the only other unknown parameter in the definition of ionization 
parameter £ = L lml /nr 2 (where L; on is the ionizing luminosity in 
the 1-1000 Ryd band). Searching for absorption line variability 
on short timescales and thus constraining the location of the WA 
is the main motivation for conducting monitoring campaigns of 
AGN (e.g., Mrk 509, see Kaastra et al. 2012). The knowledge 
of the location is crucial for estimating the mass outflow rates 
and kinetic luminosities associated to these absorbers (e.g., 
Crenshaw & Kraemer 2012) and thus, to evaluate their potential 
impact on the host galaxy environment. 

With this aim, during the summer of 2013 and the winter 
of 2013-14 we performed a large multiwavelength monitoring 
campaign on the bright Seyfert 1 NGC 5548. The overview of 
the campaign in presented in Mehdipour et al. (2015, hereafter 
Paper I). NGC 5548 is a prototypical Seyfert 1 galaxy, that 
has been studied for decades from optical (e.g., Peterson & 
Wandel 1999) to X-ray wavelengths (e.g., Nandra et al. 1993; 
Iwasawa et al. 1999). From a dynamical modeling of the BLR 
(Pancoast et al. 2014), it is inferred that this source is observed 
at an inclination angle of ~ 30° and hosts a supermassive black 
hole (SMBH) of ~ 3 x 10 7 M 0 in its center. Previously, high 
resolution UV (Crenshaw & Kraemer 1999; Crenshaw et al. 
2003a) and X-ray (Kaastra et al. 2000, 2002; Steenbrugge et al. 
2003, 2005) spectra have revealed several deep NAL that can be 
ascribed to a moderate velocity (v out =200-1200 km s' 1 ) ionized 
outflow. 

Unexpectedly, throughout the whole 2013-14 campaign 
NGC 5548 appeared dramatically different from the past (e.g., 
from the Chandra observation of 2002, Steenbrugge et al. 2005) 
being ~25 times less luminous in the soft X-rays. Moreover, it 
showed broad, asymmetric absorption troughs in the blue wings 
of the main UV broad emission lines (e.g., in Ly a, C iv, N v). In 
Kaastra et al. (2014, hereafter K14) we proposed that all these 
changes can be ascribed to the onset of a persistent, weakly 
ionized, fast (v~5000 km s' 1 ), wind (hereafter “the obscurer”). 
The obscurer is located within or just outside the BLR, at a 
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Table 1 . List of the XMM -Newton and Chandra datasets used for the present analysis. 


Obs. 

Date 

Satellites F 203oA “ 

Fo.3-2.OkeV " Fj.O-lO.OkeV ” 


yyyy-mm-dd 

(10~ 14 erg s -1 cnr 2 A~') 

(10~ 12 erg s -1 cm -2 ) 


XM1 

2013-06-22 

1 4 

2.07 

1.1 

13.2 

XM2 

2013-06-30 

1 4 

2.28 

3.5 

29.9 

XM3 

2013-07-07 

1 4 

2.07 

2.1 

21.2 

XM4 

2013-07-11 

1 4 

2.32 

3.6 

32.1 

XM5 

2013-07-15 

1 4 

2.57 

2.6 

26.9 

XM6 

2013-07-19 

1 4 

2.74 

2.3 

27.4 

XM7 

2013-07-21 

1 4 

2.88 

2.3 

23.5 

XM8 

2013-07-23 

1 4 

3.09 

2.2 

25.1 

XM9 

2013-07-25 

1 4 

2.97 

3.2 

29.8 

XM10 

2013-07-27 

1 4 

3.29 

3.0 

29.4 

XM11 

2013-07-29 

1 4 

3.21 

2.8 

26.4 

XM12 

2013-07-31 

1 4 

3.21 

2.2 

23.6 

CHI 

2013-09-01 

24 

4.54 

7.7 

38.7 

CH2 

2013-09-02 

24 

4.48 

4.4 

29.7 

CH3 

2013-09-10 

2 3 4 

4.89 

5.7 

36.1 

XM13 

2013-12-20 

1 4 

3.27 

2.1 

21.8 

XMM 

2014-02-05 

1 4 

3.49 

4.1 

24.4 


Notes. Instruments are labeled with numbers, as followed: 1: XMM -Newton, 2: Chandra , 3: NUSTAR, and 4: Swift. (<I) UV flux measured by 
Swift-UVOT, corrected for reddening . (b> Observed flux in the quoted bands, as derived from our best fit model. 


distance of a few light days from the SMBH, and possibly has 
been launched from the accretion disk. It blocks ~90% of the 
X-ray flux along our line of sight, thereby lowering the ionizing 
luminosity received by the WA. Indeed, in this obscured 
condition, the historical NGC 5548 WA is still present, but with 
a lower ionization. In the X-rays it is consistent with being ~3 
times less ionized than what was observed in 2002 (K14), and, 
at the same time, in the UV it shows new lower-ionization NAL 
(from e.g. C n and C m, Arav et al. 2014, in press). 

In this paper we use all the high resolution X-ray spectra 
collected during our campaign to assess what drives the spectral 
changes of NGC 5548 on timescale as short as few days. This 
is the typical time separation between the X-ray observations 
of the campaign. These spectra are suitable to investigate the 
absorption variability, because they cover the energy band where 
the main ionized and neutral absorption features of e.g. oxygen 
and iron fall. During the campaign the source was always weak 
in the soft X-rays, except for a sudden brightening in September 
2013 (Fig. 1). On that occasion, we triggered a Chandra -LETGS 
observation. In the following we investigate also the possible 
causes and consequences of this sudden brightening. 

The paper is organized as follows: in Sect. 2 we briefly present 
the datasets that we use in this analysis, and in Sect. 3 we 
describe the template spectral model that we apply to all the 
datasets in Sect. 4. Finally in Sect. 5 we discuss our results and 
in Sect. 6 we outline the conclusions. 

The C-statistic (Cash 1979) is used throughout the paper, and 
errors are quoted at 68% confidence level (AC = 1.0). In all 
the spectral models presented in the following, we use the 
Galactic hydrogen column density from Wakker et al. (2011, 
Ah =1.45 X 10 20 cm -2 ). The cosmological redshift that we 
adopted for NGC 5548 is 0.017175 (de Vaucouleurs et al. 
1991).The cosmological parameters are set to: Hq-70 km s 1 
Mpc -1 , D m =0.3 and Qa-0-7. 

2 . The data 

XMM-Newton observed NGC 5548 between June 2013 and 
February 2014 using both the EPIC cameras (Turner et al. 2001; 


Striider et al. 2001) and the Reflection Grating Spectrometer 
(RGS, den Herder et al. 2001). The core of the campaign 
consisted of 12, ~50 ks long, XMM -Newton observations that 
were taken every ~ 2-8 days in June and July 2013. After 
these, two other observations were acquired in December 2013 
and February 2014, providing 14 XMM -Newton datasets in to¬ 
tal. The details of all the observations, that hereafter we label 
in chronological order as XM1-XM14, are given in Table 1. 
During the entire campaign, the source was also monitored daily 
by Swift (Gehrels et al. 2004), both in the X-rays with the X Ray 
Telescope (XRT, Burrows et al. 2005) and in the UV with the UV 
Optical Telescope (UVOT, Roming et al. 2005). Swift-UVOT 
flux measurements in the UVW2 (d=2030 A) filter, corrected 
for reddening and for the host galaxy contribution as explained 
in Paper I, are also used in the present analysis. 

In September 2013 we triggered a Chandra observation because 
Swift observed a sudden brightening (Fig. 1). The observa¬ 
tion was performed with the Low Energy Transmission Grating 
Spectrometer (LETGS, Brinkman et al. 2000) in combination 
with the High Resolution Camera (HRC-S). The observing time 
was split in three observations (Obs. CH1-CH3) of 30, 67 and 
123 ks respectively. The first two were taken on September 1st 
and 2nd while the third and longest one was taken a week later on 
September 10th. Simultaneously with this observation a higher 
energy spectrum was acquired with the Nuclear Spectroscopic 
Telescope Array (NuSTAR, Harrison et al. 2013) satellite. In the 
occasion of this flaring event, no XMM -Newton observation was 
available. 

A detailed description of the data reduction procedure for all the 
instruments is given in Paper I. In the present analysis, we fit 
the high resolution RGS and LETGS spectra. The simultaneous 
EPIC-pn and NuSTAR spectra at higher, less absorbed, energies 
provide the continuum baseline over which the absorption is su¬ 
perposed. We defer the reader to Ursini et al. (2015, in press), 
hereafter U15, for a detailed modeling of the continuum at high 
energies. 

As in K14, we fitted the EPIC-pn in the 1.03-10.0 A (=1.24-12 
keV) band together with the RGS in the 5.68-38.23 A (=0.32- 
2.2 keV) band. Thus, the two instruments overlap in a small band 
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allowing to check for possible intercalibration mismatches. In 
all our fits the intercalibration factor used in K14 (~ 1.027) was 
adequate. Due to an incomplete correction for the gain of the 
EPIC-pn (not corrected in the SAS vl3 we used, see Paper I), 
the Fe Kar line appears blueshifted (see Cappi et al. in prep, for a 
detailed discussion), which we correct for with an artificial red- 
shift for these spectra. However, this solution leads to a poor fit 
near the energy of the gold M-edge of the telescope mirror. For 
this reason, we omitted the interval 5.0-6.2 A (^2.0-2.5 keV) 
from all our fits. 

We fitted the Chandra -LETGS spectra between 2 and 40 A 
(^0.3-6.2 keV). We used the NuSTAR spectrum simultaneous 
to CH3 in the 0.2-2.5 A (-5-60 keV) band. Since Chandra- 
LETGS and NuSTAR were consistent in the overlapping band 
we did not apply any intercalibration correction in the joint fit. 


3. The template model 

In order to assess the variability of NGC 5548, we used the same 
template model to fit all the spectra. This template is close to the 
model adopted in K14, differing from it only in the modeling 
of the “soft-excess” (Arnaud et al. 1985) component. We per¬ 
formed all the spectral analysis with the latest version of SPEX 
(v. 2.05, Kaastra et al. 1996). In our template model we consid¬ 
ered the cosmological redshift and the Galactic absorption. For 
the latter, we use a collisionally ionized plasma model (HOT), 
with a nominal temperature of 0.5 eV for a neutral gas. 

Our continuum model includes a primary and reflected power- 
law, on top of which lies a soft-excess component. In 
NGC 5548 the reflection is consistent to be constant (see U15) 
thus producing a steady narrow Fe KoTine (see Cappi et al., in 
prep.). In the present analysis we used the reflection parameters 
obtained in U15 (F = 1.9 and E cut =300 keV for the photon in¬ 
dex and the high energy cutoff of the primary power law) in the 
SPEX reflection model REFL, which includes both a Compton 
reflected continuum (Magdziarz & Zdziarski 1995) and the Fe 
KoTine (Zycki & Czerny 1994). To adjust the fit of the Fe 
Kg' line, we let the normalization of REFL free to vary within 
the errors of the U15 model. When using EPIC-pn data, we ap¬ 
plied to REFL an artificial blueshift (z=-1.25xlCT 2 ) to correct 
for an apparent centroid shift of the Fe Kg' line (see also Sect. 
2 ). 

The X-ray obscuration affects the spectrum mainly below 2.0 
keV and makes therefore the detection of the soft-excess from 
X-ray spectra elusive. Nevertheless, this component contributes 
to the continuum in the band where most of the absorption is 
seen, thus its modeling is critical for evaluating the absorption 
variability. In Paper I we show that the soft-excess component 
is likely to be the tail of a component extending from UV to 
soft X-rays, which is produced by Compton up-scattering of 
the disk photons in a warm, optically thick plasma. Therefore, 
in this framework, the UV emission (which is not affected by 
the obscuration) is a proxy for the X-ray soft-excess. We used 
COMT (based on Titarchuk 1994) to model the soft-excess com¬ 
ponent and we used the UV flux (F 203 0 a) listed in Table 1 for 
each observation to set the normalization. The 0.3-2.0 keV lu¬ 
minosity of the soft-excess is given by: (Ly^^^/lO^W) = 
2.093 + 2.893 X (,F 20 3 0 ^/10~ 34 Jy), valid for the range of ob¬ 
served UV fluxes used here (Mehdipour, private communica¬ 
tion). Hence, consistently with the long-term variability analysis 
in Mehdipour et al. (in prep.), we kept the COMT component 
constant in shape in all the fits, with its 0.3-2.0 keV flux vary¬ 
ing according to the UV flux. We use the mean values given in 


Paper I for the other COMT parameters, namely the Wien tem¬ 
perature of the incoming photons ( Tq = 80 eV), the temperature 
(T\ = 0.17 keV), and the optical depth (r = 21) of the plasma. 

A detailed modeling of the X-ray emission features of 
NGC 5548 will be presented in Whewell et al, in preparation. 
The X-ray narrow emission lines are consistent with being con¬ 
stant during the XMM-Newton campaign and thus could be kept 
frozen in all our fits. For the narrow lines we the values obtained 
in the K14 fit, while for the broad emission features we used the 
fluxes given in Steenbrugge et al. (2005). 

In K14 we found that the obscurer causing the persistent flux 
dimming of NGC 5548 comprises two ionization phases (here¬ 
after labeled as “warm” and “cold”). The warm phase is mildly 
ionized (log J = — 1.25) and has a larger covering fraction (e.g., 
Cv, warm ~ 86%, Cv.coid ~ 30%, K14) and a lower column 
density (e.g., V H ,warm ~ 10 22 cnT 2 , V H , C oid ~ 10 23 cm' 2 , 
K14) than the cold phase which is consistent with being neu¬ 
tral. In SPEX, we used a photoionized absorber model (XABS) 
for both the obscurer components and we adopted for them the 
same kinematics (outflow and broadening velocity) used in K14. 
The ionization balance for both the obscurer components is com¬ 
puted using the intrinsic, unabsorbed spectral energy distribution 
(SED). 

Finally, we included six more XABS components in the tem¬ 
plate model to account for the lower ionized counterpart of the 
historical NGC 5548 warm absorber. The ionization balance for 
all the WA components is computed using a SED filtered by the 
obscurer. We label the WA components with capital letters, from 
A to F as the ionization increases. As in K14, we assumed that 
the WA varies only in ionization when the ionizing SED illu¬ 
minating it changes. Indeed, from UV spectra we know that the 
kinematics of the WA components has not changed over at least 
16 years, and Ebrero et al. (in prep.) shows no historical evidence 
for total hydrogen column density variations. We kept therefore 
the kinematics and the equivalent hydrogen column density of 
all the WA components frozen to the parameters obtained from 
the unobscured 2002 spectrum (see the updated fit in Ebrero et 
al., in prep.). 

In the following we apply this template model to all the XMM- 
Newton and Chandra datasets of our monitoring campaign. In 
all the fits, the free parameters are the slope and the normal¬ 
ization of the primary power law, the column density, and the 
covering fraction of the two obscurer phases. Besides these, in 
the fits of those datasets where either the continuum or the ob¬ 
scurer parameters change significantly, we also allow the ioniza¬ 
tion parameters of the warm absorber to vary. The variability in 
flux of the soft excess component was determined directly from 
the Swift UV flux, as explained above. 

4. Fitting the variability 

4.1. The core of the campaign: from June to August 2013 

The XMM-Newton monitoring campaign began on June 22nd, 
catching NGC 5548 at the lowest flux level observed in any of 
the high resolution spectra of the campaign CF 0 . 3 - 2 . 0 kev ~ 1.1 X 
ICC 12 erg s 1 cm -2 , for Obs. XM1, see Table 1). Throughout all 
the summer of 2013, (Obs. XM2-XM12) the 0.3-2.0 keV flux 
remained quite stable (within a factor ~ 1.3) around an aver¬ 
age level of ~ 2.7 x 10~ 12 erg s _1 cm' 2 . The stacked XMM- 
Newton spectrum of Obs. XM1-XM12 is published in K14. In 
that paper we determine consistently both the average obscurer 
parameters and the average obscured SED illuminating the WA. 
Accordingly, we compute the ionization balance and the ioniza- 
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Table 2. Best fit parameters and errors, for the individual XMM -Newton and Chandra observations. In the last row, we list also the parameters 
derived in K14 for the average XM1-12 spectrum for comparison. 


Obscurer components " WA components b 


Obs 

r‘ 

M c 

1 v H, warm 

C d 

V, warm 

Nh.coIJ ‘ 

C d 

cold 





£e * 


C/Expected 



(10 22 cnf 2 ) 


(10 22 cm- 2 ) 




(log erg 

cm s 1 

) 



XM1 

1.5 ** 

1.7 ±0.1 

0.86 ±0.01 

12 ± 1 

0.47 ± 0.02 







384/334 

XM2 

1.58 ±0.03 

1.04 ±0.04 

0.927 ± 0.005 

10 ±2 

0.170 ±0.03 







422/345 

XM3 

1.5 ** 

1.42 ±0.06 

0.909 ± 0.005 

11 ±3 

0.18 ±0.02 







379/337 

XM4 

1.60 ±0.03 

1.17 ±0.04 

0.922 ± 0.004 

14 ± 2 

0.19 ±0.03 







406/346 

XM5 

1.57 ±0.03 

1.44 ±0.05 

0.915 ± 0.004 

8 ± 1 

0.24 ± 0.03 







435/335 

XM6 

1.5 ** 

1.4 ±0.1 

0.916 ±0.005 

5 ± 1 

0.30 ±0.08 

0.33* 

1.06* 

1.70* 

1.91* 

2.48* 

2.67* 

415/341 

XM7 

1.5 ** 

1.40 ±0.05 

0.909 ± 0.005 

12 ±2 

0.27 ±0.03 ‘ 

394/340 

XM8 

1.53 ±0.03 

1.62 ±0.07 

0.88 ± 0.05 

9 ± 1 

0.35 ± 0.03 







389/339 

XM9 

1.64 ±0.03 

1.58 ± 0.05 

0.896 ± 0.004 

11 ±2 

0.24 ± 0.03 







390/340 

XM10 

1.57 ±0.03 

1.43 ± 0.05 

0.911 ±0.004 

9 ± 1 

0.24 ± 0.03 







434/344 

XM11 

1.56 ±0.03 

1.46 ±0.05 

0.902 ± 0.004 

10 ±2 

0.23 ± 0.03 







399/340 

XM12 

1.5 ** 

1.43 ± 0.04 

0.922 ± 0.005 

12 ±2 

0.27 ± 0.03 







447/338 

CHI 

1.83* 

1.7 ±0.2 

0.74 ± 0.03 

17* 

< 0.23 

0.38* 

1.11 * 

1.75 * 

1.96* 

2.53* 

2.72 * 

149/152 

CH2+3 

1.83 ±0.02 

1.49 ±0.08 

0.79 ± 0.07 

17 ±3 

0.35 ± 0.03 

0.38 

1.11 

1.75 

1.96 

2.53 

2.72 

530/311 

XM13 

1.5 ** 

1.7 ±0.1 

0.87 ± 0.02 

9 ± 2 

0.33 ± 0.03 

0.33 * 

1.06* 

1.70* 

1.91 * 

2.48 * 

2.67* 

415/340 

XMM 

1.59 ±0.01 

1.22 ±0.03 

0.917 ±0.003 

10* 

< 0.003 

0.39 

1.12 

1.76 

1.97 

2.54 

2.73 

485/342 

K14 *** 

1.566 ±0.009 

1.21 ± 0.03 

0.86 ±0.02 

9.6 ±0.5 

0.30 ±0.10 

0.33 

1.06 

1.70 

1.91 

2.48 

2.67 



Notes. (a> For the warm obscurer: log£ = -1.25. For the cold obscurer: log£ = -4 (A) For Obs. CH2+3. XMM, and K14 the ionization parameters 
given are the output of the iterative fitting procedure described in Sect. 4.2. (c) Photon index of the primary continuum. <d> Column density of the 
obscurer components. w Covering fraction of the obscurer components.Ionization parameters of the warm absorber components. (? ' C-statistics 
of the final best fit model. w Frozen parameters. ( ** ) Lower limit of the fitting range. ( *** ) Best fit parameters derived in K14 for the co-added 
XM1-12 spectrum. 


tion parameters for all the WA components. 

We fitted Obs. XM1-XM12 assuming the average WA deter¬ 
mined in K14. Hence, only the continuum and the two obscurer 
phases were left free in the fits. Furthermore, in the fits we set 
a lower limit for the continuum slope (T > 1.5) to be consistent 
with what was observed during the campaign at higher energies 
and with the analysis repeated by Cappi et al, in preparation. 
Indeed, when fitting such absorbed spectra below ~10 keV, a 
flatter intrinsic continuum becomes degenerate with a stronger 
obscuration. The assumption on the slope aimed at minimizing 
this ambiguity and, however, did not affect the final results (as 
we checked a posteriori, see below). 

As a first step we fit the spectra allowing for both the obscurer 
phases just one of the parameters (column density, and cover¬ 
ing factor) free. When only the covering fractions are allowed 
to vary the fit displays strong residuals, for instance the negative 
residuals just below ~10 A (See Fig. 2, first panel). In the obser¬ 
vation where they are more evident (Obs. XM5) the C-statistic 
is 600 for an expected value of 340. When we fixed the covering 
factors but instead allowed the column densities to vary, more 
systematic residuals between 10 and 20 A are apparent (Fig. 2, 
second panel). Therefore, we conclude that the two phases of the 
obscurer have to be variable both in column density and the cov¬ 
ering fraction to adapt the template model to all the individual 
XMM -Newton observations. 

We tested however the possibility that the obscurer varies in ion¬ 
ization rather than in covering fraction. For most of the datasets, 
a statistically acceptable fit can be achieved keeping the covering 
fraction of both the obscurer phases frozen to the average val¬ 
ues derived in K14 and allowing the ionization parameter of the 
warm obscurer to vary instead. However, with these constraints 
the best fit prefers an almost neutral obscurer (e.g. log 4' < -3.5), 
which would be too lowly ionized to produce, for instance, the 
broad C iv absorption lines that are seen in the UV. Hence, we 


discarded these fits. 

The final best fit parameters are listed in Table 2. In a couple 
of cases, the fit stopped at the lower limit we had imposed for 
F. We checked how much further the fit of these datasets could 
be improved allowing an even flatter continuum. In all cases, re¬ 
leasing the spectral index resulted in a negligible improvement 
of the fit (e.g., for Obs. XM1, AC = -1 for T= 1.46). 

The RGS spectra of the individual observations are rather noisy 
and the residuals do not show any hint of unaccounted WA fea¬ 
tures. However, as a final test, we checked how sensitive are 
the best fits to possible variations in the WA ionization. We 
attempted to refit each observation setting the WA ionization 
parameters according to the variation from the average of the 
continuum normalization. All the fits were insensitive to this 
variation, (e.g., AC ~ -1) with the free parameters remaining 
the same within the errors. Therefore, we concluded a poste¬ 
riori that assuming a constant WA in the core of the XMM- 
Newton campaign was reasonable. We show in Fig. 2 (third and 
fourth panel) an example of best fit (Obs. XM5). 

4.2. The flare of September 2013 as seen with 
Chandra -LETGS and NuSTAR 

In September 2013 we triggered a series of 3 Chandra -LETGS 
observations because it seemed that NGC 5548 was recovering 
from the obscuration. Indeed, in a few days, the X-ray flux 
in both the Swift-XRT bands rose above the level measured 
at unobscured epochs and remained steady for about a week. 
This brightening was however a short-lived flare, and after a 
few days the source fell again to the typical low flux level of 
the XMM -Newton campaign (Fig. 1). Our triggered Chandra- 
LETGS observations missed the peak of the flare. The first two 
observations were taken during its declining tail, while a week 
later the third one caught a smaller rebrightening. 
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Wavelength (A) 


Fig. 2. Examples of fits the XMM -Newton spectrum of Obs. XM5. 
From the top to the bottom panel: fit residuals when only the obscur- 
ers covering fractions are allowed to vary from the values obtained for 
the average spectrum; fit residuals when only the obscurers column den¬ 
sities are allowed to vary from the values obtained for the average spec¬ 
trum; fit residuals of the best fit model where both column densities and 
covering fractions are permitted to vary freely; spectrum of observation 
XM5. The solid line represents the best fit model. The data have been 
rebinned for clarity. 


To understand if and how the absorbing components had 
responded to these continuum changes we needed to use in the 
photoionization modeling of the obscurer a SED representative 
of NGC 5548 during the flare (Fig. 3, solid line). To construct 
it, we used the Comptonization model of Paper I which extends 
from the UV to the soft X-rays. At the same time, NuSTAR 
provided the continuum slope at high energies. In the UV we 
took the model values corresponding to A — 2987 A, A — 1493 
A A = 909 A. We then interpolated 20 data points in the model 
between 0.03 and 100 keV. Finally, we cut off the SED at low 
energies (below 0.01 Ryd). 

We derived consistently the obscurer parameters and the WA 
ionization parameters using the same iterative method of K14. 
At each iteration of this fitting routine new obscurer parameters 
are fitted. Next, the new obscured continuum is used as the ioniz¬ 
ing SED in the photoionization modeling of the six WA compo¬ 
nents. The new ionization parameters are assigned by rescaling 


£ 

o 


O 

O 



E (eV) 


Fig. 3. The unobscured (solid line) and obscured (dashed line) spec¬ 
tral energy distributions for NGC 5548 during the September 2013 
flare, which were used for the photoionization modeling of the ab¬ 
sorbers. The 3 data points in the UV (filled circles) are taken from the 
Comptonization model of Paper I. The data points from the EUV to hard 
X-rays (open squares) are also interpolated from the X-ray continuum 
model. See Sect. 4.1 for details. 


those observed in the unobscured spectrum of 2002 to the level 
of the current obscured continuum. Explicitly, at the N,/, itera¬ 
tion 2002 ~ Tion,N/Tion,2002- Finally, the ionization balances 
for the new ionization parameters are recomputed before mov¬ 
ing to the next iteration. The final outputs of this procedure are 
the obscurer parameters, the obscured SED illuminating the WA, 
and the rescaled WA ionization parameters. 

At first, we dealt with the third and longest Chandra -LETGS 
spectrum (Obs. CH3) and we fitted it jointly with the simultane¬ 
ous NuSTAR spectrum. We started with a few iterations where 
only the continuum was allowed to vary. However, since in this 
way the fit was left with large residuals (C/Expected C=846/307, 
Fig. 4, first panel) we released first the covering fractions (Fig. 
4, second panel) and in turn the column densities of both the ob¬ 
scurer phases (Fig. 4, third panel). Once we achieved the best 
fit, we applied it to the other two Chandra-LETGS spectra, for 
comparison. For obs. CHI, the fit tends to steepen the contin¬ 
uum up to T ~ 2.2. In contrast, we could easily fit Obs. CH2 just 
by renormalizing the model. Therefore we decided that we can 
stack observations CH2 and CH3, and thus increase the signal- 
to-noise ratio of the spectrum. 

We fitted this stacked spectrum (hereafter labeled as Obs. 
CH2+3) together with NuSTAR using the iterative procedure 
just described. The final best fit model is shown in Fig. 4 (fourth 
panel) and the final obscured SED produced in the iterative fit¬ 
ting is plotted in Fig. 3, dashed line. Using the above best fit, 
we tested whether the difference in photon index for observation 
CHI could be due to changing properties of the obscurer. If we 
keep the continuum shape frozen, the best fit prefers zero cover¬ 
ing fraction for the cold obscurer. For this dataset, we favor this 
solution because a large variation of the continuum slope would 
be inconsistent with what NuSTAR and INTEGRAF have shown 
for the whole campaign (see U15). 

All the best fit parameters for the Chandra observations are 
shown in Table 2. Compared to the spectra of the core of the 
campaign, the Chandra spectra require both a steeper contin¬ 
uum (T ~ 1.8) and a lower covering fraction of the warm ob- 
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Fig. 4. Examples of fit for the Chandra-LETGS plus NuSTAR spec¬ 
trum of NGC 5548 during the September 2013 flare. From the top to 
the bottom panel: fit residuals when only the continuum is allowed to 
vary from the values obtained for the average spectrum; fit residuals 
when only the obscurers covering fraction is allowed to vary from the 
values obtained for the average spectrum; fit residuals for the final best 
fit model. The solid line represent the best fit model. The data have been 
rebinned for clarity. 


scurer ( Cv.warm ~ 0.8). However, in principle, it is possible that 
a change in the ionization state of the obscurer mimics a drop 
in the covering fraction. The data quality is not sufficient to fit 
the ionization parameter, therefore we refitted the spectrum with 
the ionization parameter of the warm obscurer expected if it re¬ 
sponds immediately to flux changes. We used the UVW2 flux 
to calculate the expected increase in ionization parameter. In the 
fit we used the covering fractions as given in K14. With these 
constraints the resulting fit is statistically worse (C/Expected 
C=647/311) and shows larger positive residuals in the Chandra- 
LETGS band. Leaving the ionization parameter of the warm ob¬ 
scurer free, a better fit is derived, but the obtained ionization pa¬ 
rameter (log £ ~ 0.9) is unrealistically high (~100 times than the 
average value) considering the increase by only a factor 2 in the 
UVW2 flux during the flare. Therefore, we can exclude that a 
change in the ionization of the obscurer is the dominant cause of 
the observed flare. 

During the XMM -Newton campaign, the discrete features of 
the WA are always blended with the obscured continuum and 



Fig. 5. The CAa/idra-LETGS spectrum of NGC 5548 in the 15-24 A 
wavelength region. The solid line represents our best-fit model. The 
most prominent emission and absorption features are labeled. The spec¬ 
trum shows some WA signatures (Fe-UTA, O iv-O v) that were not de¬ 
tected during the XMM -Newton campaign. 

Table 3. List of the O iv and O v lines predicted by our WA model that 
contributes to the features visible in Fig. 5. 


Ion 

WA component 

Outflow velocity T 0 b s “ 

km s ^ 1 A 

T 6 


C 

1148 

20.23 

27 


B 

547 

20.27 

7 


D 

254 

20.29 

3 


E 

792 

20.26 

1 


B 

547 

21.04 

3 


C 

1148 

21.00 

2 


C 

1148 

22.67 

143 

Ov 

B 

547 

22.71 

34 


D 

254 

22.73 

16 


C 

1148 

23.04 

10 


B 

547 

22.08 

5 


Notes. (a) Predicted wavelength, considering the cosmological redshift 
and the blueshift due to the outflow. (b) Line optical depth. 


not detectable. During the Chandra observation the source was 
about twice as bright, and some WA signatures are visible in 
the stacked CH2+3 spectrum (Fig. 5). These features are consis¬ 
tent with the WA model computed via our iterative procedure. 
The broad trough at ~ 16 A is a blend of unresolved transition 
array (UTA) from several ionized iron species Between 20 and 
24 A some Oiv-Ov absorption lines are present. In table 3, we 
list which of the O iv and O v lines predicted by our WA model 
contribute to each feature. The WA comprises 6 ionization com¬ 
ponents that could be distinguished thanks to the excellent data 
quality of the 2002 spectrum. Here, the lower statistics does not 
allow to overcome the blending among the components. 

4.3. The end of the campaign: the observations of December 
2013 and February 2014 

In the last two XMM-Newton observations of the campaign, 
NGC 5548 was again at the same flux level of summer 2013. 
Therefore, at first we attempted to fit them using again the same 
WA of the average spectrum. The parameters of the contin¬ 
uum and of the two phases of the obscurer were free. This at- 
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Fig. 6. XMM-Newton spectrum and model residuals for Obs. 14. The 
solid line represent our best fit model. The data have been rebinned for 
clarity. 


tempt resulted in a satisfactory fit for XM13, while for Obs. 
XM14, we obtained a high column density for the cold obscurer 
( ^h, cold ~ 7x 10 23 cm~ 2 ) largely inconsistent with what was ob¬ 
served throughout the campaign. Thus, we attempted to refit this 
dataset freezing the cold column density to the average value 
measured in the previous observations. This resulted in negli¬ 
gible covering fraction for the cold obscurer ( Cy.coid ~ 10 2 ) 
and in a flatter continuum (T ~ 1.6 instead of ~1.7). This lower 
value of T is also similar to the other values measured at this lu¬ 
minosity (Fig. 8). For these reasons, we considered this solution 
more physically plausible. We note that for this dataset a solu¬ 
tion with a negligible covering fraction of the cold obscurer and 
a flatter continuum is found also in U15 when including in the 
fit INTEGRAL data at higher energies. 

Considering the large change in the obscurer, also for this obser¬ 
vation we used our iterative method to determine the ionization 
of the WA (for details see Sect. 4.1). We built the SED for Obs. 
XM14 interpolating the Comptonization model of Paper I re¬ 
peating all the steps already described in Sect. 4.1. In this case, 
the slope of the continuum at X-ray energies is provided by the 
fit of the EPIC-pn data. Starting from the solution just described, 
the fit converged in a few iterations. We also made some itera¬ 
tions where the ionization parameter of the warm obscurer was 
left free. However, since it was not possible to constrain it, we 
finally kept it frozen to the value found in K14 (see Sect. 3). We 
adopt the final result of the iterative routine as our best-fit model. 
The parameters resulting from this exercise are outlined in Table 
2 . 

We show the final best-fit model for Obs. XM14 in Fig. 6. The 
final WA model that we found is only slightly different from the 
average one. The residuals of the best-fit model display some 
possible features at ~ 10.3 A and ~12.3 A which are close to the 
expected wavelength of the Neix-K edge (10.54 A) and NeX- 
Ly a (12.34 A) respectively. We checked if adding additional ab¬ 
sorbing column density in Ne ix-Ne x could improve the fit. For 
this, we used the SLAB model in SPEX, that allows us to fit a 
single ionic column density regardless of the ionization balance 
with the other ions. We found however that this additional com¬ 
ponent is not required by the fit (AC = -1). 


5. Discussion 

5 . 1. The short-term variability 

During our extensive monitoring campaign in 2013 and early 
2014 NGC 5548 was always obscured. In this analysis we have 
applied the model developed in K14 for the average spectrum 
of the core of the campaign to all the individual observations, 
with the aim of understanding how the source varies. When 
both the intrinsic continuum and the obscurer are allowed to 
vary, the model is able to explain the variability on the 2 days-8 
months timescale sampled in the monitoring campaign. The 
obscuring material that is causing the soft X-ray flux depression 
of NGC 5548 varies along our line of sight, both in column 
density and in covering fraction. The scenario proposed by K14, 
that the source is obscured by a patchy wind, is consistent with 
our variability findings. In this framework, the variability of the 
obscuration may well be due to several reasons, e.g. motion 
across the line of sight and changing ionization with continuum 
variability. 

The best fit parameters for the continuum and the two phases of 
the obscurer are displayed in Fig.7 as a function of time. During 
the core of the XMM-Newton campaign (Obs. XM1-XM12) 
the source was steadily obscured and the variability in flux 
was relatively small (~ 27%). The only clear outlier with a 
flux significantly different from the average is Obs. XML The 
values we found for the warm obscurer parameters (Fig. 7 4th 
and 6th panel) deviate from those found in K14 the co-added 
XM1-XM12 spectrum. This is due to our different modeling of 
the soft-excess, that dominates the continuum in the band where 
the absorption from this component is more effective. In our 
modeling the Comptonized soft-excess, whose normalization 
is set by the UV flux measured by Swift for each observation, 
is always more luminous (To^% keV = [0.8 - 1.4] x 10 43 
erg s 1 ) than the phenomenological blackbody fitted in K14. 
Throughout the core of the campaign, the intrinsic continuum 
is fairly constant in shape (with a standard deviation of the 
spectral index <x r ~ 3%) and slightly variable in normalization 
(cr uvw 2 ~ 17%, cr Norm ~ 24% for the soft-excess and the 
power law component, respectively). For the obscurer, the 
cold component was the most variable (cr co[d ~ 31% and 
<t l y H o]d ~ 23%). In contrast, for the warm component the 
covering fraction is stable (cr r ~ 2%) and the column 

° V, warm 

density shows rather small variability (cr ~ 13%). The 

large deviation from the average of the cold covering fraction 
suggests that the obscurer inhomogeneity, that is possibly 
dominated by the cold phase, may have caused most of the 
variability observed during this phase of the campaign. 

In this paper we have also presented the Chandra- LETGS 
datasets that were acquired in September 2013 when 
NGC 5548 underwent a two-week brightening. With re¬ 
spect to the core of the campaign changes in both the continuum 
and the obscurer are required to fit these spectra. At the time 
of the Chandra observation, the UV flux measured by Swift, 
which in our interpretation is a tracer for the soft X-ray excess, 
was the highest of the whole campaign. At the same time, the 
continuum at hard X-ray energies increased in flux and became 
steeper. For both observations the warm obscurer component has 
a significantly lower covering fraction. In the first observation 
also the covering fraction of the cold component is lower. As 
pointed out in Sect. 4.2, a variation in the obscurer ionization 
alone is insufficient to explain to observed variation in spectral 
shape. 
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Time (MJD-56465) 

Fig. 7. From the top to the bottom panel: the observed flux in the 0.3-2.0 keV band, the photon index of the primary continuum, the column 
densities of the cold and the warm obscurer, and the covering fractions of the cold and warm obscurer are shown as a function of time. The flux 
is plotted in units of ICE 11 erg s _1 citT 2 . The column densities are plotted in units of 10 22 cirT 2 . The Modified Julian Day (MJD) correspondent 
to each observation is labeled on the horizontal axes, which has been shrunk for display purpose. Black diamonds and red asterisks identify 
parameters measured with XMM -Newton and Chandra-LETGS, respectively. Error bars, when larger than the size of the plotting symbol, are also 
shown. Upper limits are plotted as an arrow. Crosses represent values that were kept frozen in the fits. In each panel, the dotted horizontal line 
indicates the parameter value derived in K14 for the co-added XM1-12 spectrum. 
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The decrease in covering fraction, in principle, can be due 
either to a local thin patch of the obscurer passing in our line of 
sight at the moment of the Chandra or to a geometrical change 
of the UV/soft X-ray continuum source behind the obscurer. 
The former possibility, although it cannot be excluded, would 
however require that the intrinsic continuum and the obscurer 
change properties in synchrony, which seems ad hoc. In the 
Comptonization model of Petrucci et al. (2013), the UV/soft 
X-ray spectrum is supposed to be produced via Comptonization 
of the UV disk photons in a “warm” (T ~ 1 keV) and moder¬ 
ately thick (t ~ 10 - 20) corona. A “hot” corona, with higher 
temperature (~ 100 keV) and smaller optical depth (~ 1), will 
in turn Compton upscatter these UV/soft X-ray photons to hard 
X-ray energies. In this interpretation, an increase in physical 
size of the warm corona, while naturally augmenting both the 
UV and the soft X-ray flux, would also result in a drop of the 
observed obscurer covering fraction. Moreover, the increase in 
the UV/soft X-ray photon flux will more effectively cool the hot 
corona, hence producing a steeper hard X-ray spectrum which 
is in agreement with the observation. 

Due to increased soft X-ray flux in the Chandra observations 
some discrete WA features (Fe UTA, O iv-O v) became evident 
in the spectrum. These are the only detectable WA signatures 
in any X-ray spectrum of our campaign. These features are 
best fitted by a WA which has a significantly lower degree 
of ionization than what is observed in the unobscured 2002 
spectrum. Like K14, for all the WA components we found best 
fit ionization parameters which are 0.40 dex lower than the 
2002 values ( log ^002 = 0 .78, j 51- 2 .15, 2.36, 2.94, and 3.13). 
This means that the ionizing luminosity received by the WA 
decreased by a factor of ~ 4. Thus, our analysis confirms the 
K14 finding. The decrease in the WA ionization, that is seen 
also in the UV (Arav et al. 2014) is explained when the newly 
discovered obscurer is located between the nucleus and the 
WA. In this geometry, the obscurer shadows the central source 
and prevents most of the ionizing flux from reaching the warm 
absorber. 

The absorbers in NGC 5548 changed again in the last observa¬ 
tion of the campaign, namely Obs. XM14. We found that in this 
dataset, the covering fraction of the cold obscurer became negli¬ 
gible. At the same time, the continuum above 2.0 keV is similar 
to what is observed throughout the XMM-Newton campaign 
while the soft-excess component is only slightly higher. Thus, 
the most likely cause of the spectral changes observed in this 
dataset is again the inhomogenity of the obscuration. 

Variability in the continuum and obscurer parameters is also 
noticed in the U15 analysis of the six XMM-Newton and one 
Chandra observations that were acquired simultaneously with 
a higher energy observation. However, the parameters obtained 
fitting the EPIC-pn jointly with RGS as done here are not 
directly comparable with those obtained fitting the EPIC-pn 
jointly with NuSTAR and/or INTEGRAL, as done in U15. 
This is both because of cross-calibration issues between the 
instruments and of differences in the analysis. In particular, 
as noticed also in U15, RGS and EPIC-pn are mismatched 
in flux in the overlapping band as a function of energy (e.g., 
Detmers et al. 2009). On the other hand, NuSTAR spectra 
are systematically steeper than EPIC-pn spectra (AT ~ 0.1 
see Cappi et al., in prep). In the present analysis we consider 
also the ionized phase of the obscurer (K14), while U15 use 
two purely neutral components. This can affect the broadband 
curvature of the model. Moreover U15 has an additional degree 
of freedom in the high-energy cutoff of the continuum. For all 
these reasons, the only meaningful comparison is among the 



Fig. 8. Fitted continuum slope as a function of the observed 5.0-10.0 
keV flux for all the XMM-Newton and Chandra datasets, with errors. 
Black diamonds and red asterisk identify parameters measured with 
XMM-Newton and Chandra-LETGS, respectively. The black cross rep¬ 
resents the higher value of T that, for Obs. XM14, would be required by 
a fit including the cold obscurer (see Sect. 4.3). 


overall trend of the parameters. Even taking into account the 
differences between our analysis and the one presented in U15 
the parameters trends that we discuss below still hold. 


5.2. What drives the variability? 

To understand if there are some systematic factors driving the 
short-term variability of the source, we looked for correlations 
among the best fit parameters and the unobscured flux measured 
for the 16 datasets presented here. We used the hard X-ray flux 
in the 5.0-10.0 keV band and the UVW2 flux listed in Table 1 as 
tracers of the intrinsic continuum, as they are almost unaffected 
by the obscuration. Even considering these unabsorbed bands, 
the range of flux sampled in the monitoring campaign is nar¬ 
row (a factor of —2), with the only outliers at lower and higher 
flux being Obs. XM1 and Obs CH2+3. Therefore, to evaluate 
the reliability of any correlations we checked if it is still holding 
when removing these two data points from the computation of 
the Pearson correlation coefficient. 

In Fig. 8 we show that the best fit continuum slope steepens as 
the hard X-ray flux increases. For the complete sample the cor¬ 
relation is extremely significant. The Pearson correlation coeffi¬ 
cient is R a ii = 0.85 implying a probability p ~ 10 5 for the null 
hypothesis. When excluding Obs. XM1 and Obs. CH2+3 from 
the computation, the degree of correlation still remains signifi¬ 
cant (R X m2-xmi4=0.68, Rxm2-xmi4=1%)- This trend has been al¬ 
ready noticed in the past for NGC 5548 in Kaastra et al. (2004) 
and has been also reported in other Seyfert galaxies (e.g., MCG 
6-30-15 Shih et al. 2002), with different interpretations (see e.g., 
Ponti et al. 2006; Giacche et al. 2014). In the same Fig. 8 we 
show also that for Obs. XM14, the higher value of T that would 
be required by a fit including a thick cold obscurer (that we re¬ 
jected, see. Sect. 4.3) is inconsistent with the correlation valid 
for all the other datasets. 

In Fig. 9 we plot the parameters of both the obscurer phases as 
function of the UVW2 flux. The only parameter showing a possi¬ 
ble trend with the intrinsic continuum is the warm covering fac¬ 
tor. Namely, the drop observed during the Chandra observations 
may be the tail of a mild decreasing trend visible also for 
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the XMM -Newton data points (Fig. 9, bottom panel). Formally, 
when the two Chandra data points are included in the computa¬ 
tion a significant correlation (/? a n = -0.73, p a ii=l%) is present. 
When considering only the XMM -Newton sample, the trend is 
only qualitative (/?xmi-xmi4 = -0.03 /?xmi-xmi 4 = 25 %). In Sect. 
5.1 we have suggested that the drop in covering fraction ob¬ 
served during the September 2013 flare is due to an increase in 
the size of the soft X-ray/UV source. A clear correlation between 
the warm covering fraction and intrinsic continuum, supported 
by more numerous data points at different flux values, would 
favor the hypothesis that this is a systematic effect producing 
at least part of the observed covering fraction variability. This 
trend is not apparent for the cold covering fraction. This could 
be due to a higher degree of inhomogenity in the cold phase that 
would also explain its larger variability in covering fraction (e.g., 
it went from 0.47 in Obs XM1 to 0 in Obs. XM14). 

In conclusion, a combination of changes in the continuum and in 
the obscurer physical parameters is required to explain the short¬ 
term spectral variability of NGC 5548 during our 2013-2014 
campaign. The lack of correlation between the intrinsic contin¬ 
uum (as traced by the UVW2 flux) and obscurer parameters in¬ 
dicate that the obscurer must physically change properties inde¬ 
pendent on the source flux level. The case of the September 2013 
spectrum suggests that the soft X-ray emitting region change ge¬ 
ometry as the flux increases. This could be a systematic effect 
contributing to the overall covering fraction variability. 

6. Summary and conclusions 

During the multiwavelength monitoring campaign that we per¬ 
formed in 2013-2014 for the Seyfert 1 galaxy NGC 5548, the 
source had a soft X-ray flux well below the long-term average, 
except for a two-week long flare in September 2013. In K14, we 
have ascribed this condition to the onset of a persistent, weakly- 
ionized but high-velocity wind that blocks ~90% of the soft X- 
ray flux and lowers the ionizing luminosity received by the WA. 
Thus, in this condition, the normal WA that was previously ob¬ 
served in this source is still present, but with a lower ionization. 
In this paper we fitted all the high-resolution XMM -Newton and 
Chandra datasets that were taken during the campaign with a 
model that consistently accounts for a variable continuum, the 
newly discovered obscurer and the new ionization conditions of 
the historical WA. We found that: 

1. On the timescales sampled in the monitoring campaign (2 
days-8 months) both the intrinsic continuum and the ob¬ 
scurer are variable. The obscuring material varies both in 
column density and in covering fraction. This rapid variabil¬ 
ity is consistent with the picture of a patchy wind proposed 
by K14. 

2. The Chandra spectra that were taken just after the peak of 
the flare in September 2013 are explained by both an increase 
and a steepening of the intrinsic continuum and a drop in 
the obscurer covering fraction. The latter is likely to be due 
to a geometrical change of the soft X-ray continuum source 
behind the obscurer. 

3. The Chandra spectra of September 2013 show absorption 
from Fe-UTA, O iv and O v, consistent with belonging to the 
lower-ionized counterpart of the historical NGC 5548 warm 
absorber. These are the only individual WA features in any 
X-ray spectrum of the campaign. 

4. A positive correlation between the X-ray continuum slope 
and the observed 5.0-10.0 keV flux holds for both the XMM- 
Newton and Chandra datasets. 


5. The addition of the two Chandra points produce a formal 
anticorrelation between the warm obscurer covering fraction 
and the intrinsic continuum luminosity, as traced by the ob¬ 
served UVW2 flux. 
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values that were kept frozen in the fits. 


Kaastra, J. S.. Mewe, R., & Nieuwenhuijzen, H. 1996, in UV and X-ray 
Spectroscopy of Astrophysical and Laboratory Plasmas, ed. K. Yamashita & 
T. Watanabe, 411—414 

Kaastra, J. S., Raassen, A. J. J., Mewe, R., et al. 2004, A&A, 428, 57 
Kaastra, J. S., Steenbrugge, K. C., Raassen, A. J. J., et al. 2002, A&A, 386, 427 
Kraemer, S. B., Crenshaw, D. M., Dunn, J. P., et al. 2012, ApJ, 751, 84 
Kraemer, S. B., Crenshaw, D. M., George, I. M., Gabel, J. R., & NGC 4151 


Team. 2005, in Bulletin of the American Astronomical Society, Vol. 37, 
American Astronomical Society Meeting Abstracts, 1190 
Krongold, Y., Nicastro, F., Brickhouse, N. S., Elvis, M., & Mathur, S. 2005, ApJ, 
622, 842 

Krongold, Y., Nicastro, F., Elvis, M., et al. 2007, ApJ, 659, 1022 
Longinotti, A. L., Costantini, E., Petrucci, P. O., et al. 2010, A&A, 510, A92 
Longinotti, A. L., Krongold, Y., Kriss, G. A., et al. 2013, ApJ, 766, 104 


12 





















L. Di Gesu et al: Anatomy of an AGN in NGC 5548 


Magdziarz, P. & Zdziarski, A. A. 1995, MNRAS, 273, 837 
Marchese, E., Braito, V., Della Ceca, R.. Caccianiga, A., & Severgnini, P. 2012, 
MNRAS, 421, 1803 

Matt, G., Guainazzi, M., & Maiolino, R. 2003, MNRAS, 342, 422 
McHardy, I. M., Koerding, E., Knigge, C., Uttley, P., & Fender, R. P. 2006, 
Nature, 444, 730 

Mehdipour, M., Kaastra, J. S., Kriss, G. A., et al. 2015, ArXiv e-prints 
Miniutti, G., Sanfrutos, M„ Beuchert, T., et al. 2014, MNRAS, 437, 1776 
Nandra, K., Fabian, A. C., George, I. M., et al. 1993, MNRAS, 260, 504 
Pancoast, A., Brewer, B. J„ Treu, T„ et al. 2014, MNRAS, 445, 3073 
Peterson, B. M. & Wandel, A. 1999, ApJ, 521, L95 
Petrucci, P.-O., Paltani, S„ Malzac, J., et al. 2013, A&A, 549, A73 
Piconcelli, E., Bianchi, S., Guainazzi, M., Fiore, F., & Chiaberge, M. 2007, 
A&A, 466, 855 

Ponti, G., Miniutti, G„ Cappi, M., et al. 2006, MNRAS, 368, 903 
Ponti, G., Papadakis, I., Bianchi, S., et al. 2012, A&A, 542, A83 
Puccetti, S„ Fiore, F„ Risaliti, G„ et al. 2007, MNRAS, 377, 607 
Risaliti, G., Elvis, M„ Bianchi. S„ & Matt, G. 2010, MNRAS, 406, L20 
Risaliti, G., Elvis, M., Fabbiano, G., Baldi, A., & Zezas, A. 2005, ApJ, 623, L93 
Risaliti, G., Salvati, M„ Elvis, M„ et al. 2009, MNRAS, 393, LI 
Roming, P. W. A., Kennedy, T. E., Mason, K. O., et al. 2005, Space Sci. Rev., 
120,95 

Shih, D. C„ Iwasawa, K„ & Fabian, A. C. 2002, MNRAS, 333, 687 
Steenbrugge, K. C., Fenovcfk, M., Kaastra, J. S., Costantini, E., & Verbunt, F. 
2009, A&A, 496, 107 

Steenbrugge, K. C., Kaastra, J. S„ Crenshaw, D. M„ et al. 2005, A&A, 434, 569 
Steenbrugge, K. C., Kaastra, J. S„ de Vries, C. P., & Edelson, R. 2003, A&A, 
402,477 

Striider, L., Brief U., Dennerl, K., et al. 2001, A&A, 365. L18 
Titarchuk, L. 1994, ApJ, 434, 570 

Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27 
Turner, T. J„ Miller, L., Kraemer, S. B., Reeves, J. N., & Pounds, K. A. 2009, 
ApJ, 698, 99 

Ursini, F., Boissay, R., Petrucci, P.-O., et al. 2015, ArXiv e-prints 
Wakker, B. P., Lockman, F. J., & Brown, J. M. 2011, ApJ, 728, 159 
Walton, D. J„ Risaliti, G., Harrison, F. A„ et al. 2014, ApJ, 788, 76 
Weymann, R. J., Morris, S. L.. Foltz, C. B., & Hewett, P. C. 1991, ApJ, 373, 23 
Zycki, P. T. & Czerny, B. 1994, MNRAS, 266, 653 


13 



